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ABSTRACT 

VLA  observations  at  20  and  90  cm  of  five  distant  Galactic  H  n  regions  in  an  area  including  the  Sharpless 
S53  complex  have  been  used  to  measure  their  effective  sizes  and  optically  thick  continuum  flux  densities. 
These  data  have  been  combined  with  shorter  wavelength,  optically  thin  flux  density  measurements  and  elec¬ 
tron  temperature  estimates  to  determine  accurate  emission  measures  where  previous  results  gave  only  lower 
limits.  Our  analysis  indicates  that  existing  lower  limits  to  the  emission  measures  often  underestimate  the  true 
values  by  more  than  an  order  of  magnitude. 

Filling  factors  for  the  ionized  gas  in  these  H  n  regions  have  also  been  determined  and  values  ranging  from 
0.03  to  0.50  indicate  significant  clumping.  If  these  clumps  are  embedded  in  a  lower  density  but  more  widely 
distributed  ionized  component,  even  our  much  higher  emission  measures  are  still  lower  limits  to  the  true 
values  for  the  clumps. 

Subject  headings:  interferometry  —  nebulae:  H  n  regions  —  nebulae:  individual  (Sharpless  53) 


I.  INTRODUCTION 

Since  early  measurements  by  Mezger  and  Hoglund  (1967), 
radio  continuum  and  recombination  line  observations  have 
been  used  to  estimate  electron  temperatures  for  Galactic  H  n 
regions  (see  Wilson  1980  and  references  therein).  If  an  H  II 
region  is  of  uniform  density,  its  local  thermodynamic  equi¬ 
librium  (LTE)  electron  temperature  (TJ  is  directly  related  to 
the  ratio  of  the  measured  recombination  line  ( TL )  to  continuum 
(Tc)  brightness  temperature.  The  value  of  Te  can  then  be  com¬ 
bined  with  Tc  to  determine  the  emission  measure  (EM)  of  the 
region  if  the  source  is  resolved  by  the  telescope  beam.  If  a 
region  is  either  unresolved  or  significantly  clumped,  the  EM 
determined  in  this  way  is  only  a  lower  limit. 

Most  early  analyses  were  limited  to  assuming  LTE  and 
internal  homogeneity  to  obtain  a  single  value  of  Te  and  EM  for 
each  H  n  region.  However,  over  the  last  few  years,  efforts  have 
concentrated  on  making  detailed  models  of  H  n  regions  which 
include  gradients  in  temperature  and  density  as  well  as  depar¬ 
tures  from  LTE  (Wilson  and  Jager  1987).  These  models  have 
generally  been  applied  to  well  resolved,  nearby  H  ii  regions 
such  as  the  Orion  Nebula  where  multifrequency  observations 
are  available.  For  more  distant  H  n  regions  where  such 
detailed  observations  are  lacking,  the  assumptions  of  uni¬ 
formity  and  LTE  must  still  be  made.  Such  oversimplifications 
can  lead  to  artificially  low  EM  values. 

In  this  paper  we  use  90  cm  continuum  observations  with  the 
VLA2  of  five  distant  H II  regions  near  the  Sharpless  53  complex 
to  obtain  improved  estimates  of  their  physical  properties.  Spe¬ 

1  National  Research  Council/NRL  Cooperative  Research  Associate. 

2  The  VLA  (Very  Large  Array)  is  a  facility  of  the  National  Radio 
Astronomy  Observatory  which  is  operated  by  Associated  Universities,  Inc. 
under  contract  to  the  National  Science  Foundation. 


cifically,  we  show  that  continuum  observations  at  long  wave¬ 
lengths  where  the  H  ii  regions  are  optically  thick  can  be  used 
together  with  Te  values  and  flux  density  measurements 
obtained  at  shorter  wavelengths  where  the  regions  are  opti¬ 
cally  thin  to  determine  accurate  EM.  We  also  show  that  such 
observations  can  be  used  to  constrain  the  filling  factors  of  the 
ionized  gas  in  the  H  ii  regions. 

II.  OBSERVATIONS 

Observations  of  the  area  around  the  H  ii  complex  S53  were 
obtained  at  1480.7  MHz  and  1536.2  MHz  (220  cm,  12  antenna 
subarray)  using  50  MHz  bandwidth  and  at  327.5  MHz  and 
333.0  MHz  (290  cm,  15  antenna  subarray)  using  3.0  MHz 
bandwidth  during  1987  March  in  the  “D”  configuration 
(maximum  baseline  ~  1  km)  of  the  VLA.  The  individual  maps, 
after  calibration,  correction,  and  inspection  for  quality,  were 
added  to  produce  a  single  map  for  each  observing  band.  These 
two  “dirty”  maps  were  then  CLEANed  and  restored  with 
elliptical  Gaussians  to  an  effective  half-power  beamwidth 
(HPBW)  of  312  x  218  at  position  angle  25°  at  90  cm  and  L6 
x  0!9  at  position  angle  —21°  at  20  cm.  The  final  maps,  after 
correction  for  primary  beam  attenuation,  are  shown  in  Figures 
1  (20  cm)  and  2  (90  cm). 

The  total  integration  time  for  each  map  was  approximately 
1  hr  and  consisted  of  six  10  minute  snapshots  spread  over  an 
8  hr  period  in  order  to  maximize  the  available  uv  coverage. 
The  measured  rms  noise  level  on  the  20  cm  map  is  ~  0.005  Jy 
beam-1,  and  the  corresponding  figure  on  the  90  cm  map  is 
~0.05  Jy  beam-1. 

Calibration  at  20  cm  used  single  20  minute  observations  of 
primary  calibrators  3C  286  and  3C  48,  defined  to  have  flux 
densities  of  14.51  and  15.76  Jy,  respectively,  to  determine  the 
variable  flux  density  of  1741—038,  the  secondary  calibrator. 
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Fig.  1. — Final  contour  map  at  20  cm  of  the  region  around  the  H  n  complex  S53.  Five  distinguishable  H  n  regions  are  arbitrarily  labeled  A-E,  and  an  adjacent 
nonthermal  source  (see  text)  is  labeled  F.  The  positions  of  HI  10a  detections  (Downes  et  al.  1980)  are  marked  with  plus  signs.  Contour  levels  are  (0.5, 1, 2,  3, .. .,  10, 15, 
20, 25, . . .)  x  0.05  Jy  beam  ~ 1  with  a  peak  brightness  of  2.2  Jy  beam  ‘ 1  at  source  C.  An  ellipse  illustrating  the  synthesized  half-power  beam  width  (HPBW)  is  shown  in 
the  lower  right-hand  corner. 


RIGHT  ASCENSION 

Fig.  2. — Final  contour  map  at  90  cm  of  the  region  around  the  H  n  complex  S53.  The  distinguishable  sources  (B-F)  are  labeled  as  in  Fig.  1  with  source  A  no  longer 
detectable  at  this  frequency.  The  source  at  the  right-hand  edge  of  the  map  near  a.  =  18b21m10s,  5  =  — 13°07'  is  apparently  nonthermal  and  is  thought  to  be  unrelated 
to  the  complex.  The  contour  levels  are(l,  2,4, 6,  8, 10, 15, 20,  30, 40, . ..)  x  0.05  Jy  beam" 1  with  a  peak  brightness  of  2.0  Jy  beam"1  at  source  F.  An  ellipse  illustrating 
the  synthesized  half-power  beam  width  (HPBW)  is  shown  in  the  lower  right-hand  comer. 
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Observations  for  ~10  minutes  each  hour  of  1741—038  then 
were  used  to  fix  the  phase  corrections  and  flux  density  scale  for 
the  20  cm  data  set.  At  90  cm  no  secondary  calibrators  were 
used,  and  flux  density  and  phase  calibration  were  based  on 
observations  of  3C  286,  3C  380,  and  3C  48,  which  were  defined 
to  have  flux  densities  of  29  Jy,  45  Jy,  and  47  Jy,  respectively. 
Standard  computer  programs  were  used  to  interpolate  the  gain 
and  phase  of  the  instrument  between  calibrator  observations. 
These  procedures  are  normal  for  the  VLA  and  are  thought  to 
provide  a  phase  and  flux  density  calibration  which  is  consistent 
to  a  few  percent  with  possible  systematic  errors  of  up  to  ~  5% 
at  both  wavelengths.  For  a  more  detailed  description  of  the 
VLA  see  Napier,  Thompson,  and  Ekers  (1983). 

The  shortest  array  spacings  at  20  cm  and  90  cm  were  ~  2002 
and  ~  552,  respectively.  This  means  that  sources  larger  than 
~  3'  at  20  cm  and  ~  10'  at  90  cm  may  be  undersampled  and 
have  their  flux  densities  underestimated  by  >  ~  10%. 
However,  since  our  individual  sources  are  all  <  ~4'  in  size 
(except  source  E  which  is  ~  7'  in  one  dimension),  this  is  not  a 
problem,  especially  not  for  the  more  critical  90  cm  flux  den¬ 
sities.  Evidence  for  the  accurate  detection  of  large-scale  struc¬ 
ture  is  also  given  by  the  fortuitous  presence  of  the  extended 
SNR  G18.8  +  0.3  within  the  90  cm  field  of  view.  This  source  is 
much  larger  than  sources  A-E  and  is  comparable  in  size  (~  25') 
to  the  whole  S53  complex.  Its  radio  spectrum  is  well  deter¬ 
mined  from  centimeter  to  dekametric  wavelengths  (see  Kassim 
1987)  and  our  90  cm  integrated  flux  density  of  37  Jy  is  only 
~25%  lower  than  the  expected  50  Jy.  This  is  well  within  the 
~  50%  error  to  which  the  90  cm  flux  density  can  be  predicted 
from  the  known  spectrum. 

A  difficulty  with  90  cm  observations  at  the  VLA  is  crosstalk 
between  antennas  of  the  array.  This  is  particularly  severe  in  the 
“D”  configuration  for  low-declination  sources,  when  the 
antennas  are  pointed  nearly  towards  the  backs  of  their  neigh¬ 
bors,  and  can  lead  to  large,  low-frequency  “  ripples  ”  across  the 
maps.  Intensive  investigation  of  the  problem  and  consultation 
with  other  observers  showed  that  such  ripples  could  be  mini¬ 
mized  by  eliminating  many  of  the  short  baselines  in  the  obser¬ 
vations.  In  our  case,  this  was  found  to  mainly  affect  telescopes 


on  the  southeast  arm  of  the  VLA,  and  the  shortest  baselines  on 
this  arm  had  to  be  eliminated  from  the  data  set.  Even  though 
the  maps  still  do  not  reach  the  theoretical  random  noise  of 
~ 0.002  Jy  beam-1  expected  for  the  90  cm  observations,  no 
other  data  editing  or  correcting  procedures  could  be  found  to 
further  improve  the  results.  As  time  and  money  permit,  shield¬ 
ing  is  being  installed  on  individual  telescopes  at  the  VLA  to 
reduce  or  eliminate  this  problem.  In  the  meantime,  our  editing 
reduced  the  effects  to  manageable  levels  while  still  permitting 
the  retention  of  enough  short  90  cm  baselines  to  retain  large 
scale  structure  information  in  the  maps. 

III.  RESULTS 

The  20  cm  map  (Fig.  1)  reveals  six  discrete  sources  which 
have  been  arbitrarily  labeled  A-F  for  reference.  Plus  signs 
show  the  positions  where  the  6  cm  HI  10a  recombination  line 
emission  has  been  found  (Downes  et  al.  1980)  thus  identifying 
sources  A-E  as  H  ii  regions.  Source  F  shows  no  HI  10a  emis¬ 
sion  and  also  has  a  steep  spectrum  (see  below)  so  that  it  is 
almost  certainly  nonthermal. 

The  parameters  for  the  HI  10a  line  observations  are  listed  in 
the  upper  part  of  Table  1.  They  include  measurements  of  the 
brightness  temperature  in  the  continuum  and  in  the  line,  the 
velocity  of  the  line,  and  estimates  of  EM  and  Te  calculated  by 
assuming  the  sources  are  optically  thin  and  in  LTE.  These 
estimates  of  EM  are  lower  limits  for  sources  unresolved  by  the 
2'J  beam  of  the  Effelsberg  100  m  telescope  or  for  sources  with 
significant  clumping.  For  sources  A  and  C  the  parameters 
(except  for  the  EM)  are  also  given  in  the  lower  part  of  Table  1 
from  the  Effelsberg  H76a  recombination  line  survey  by  Wink, 
Wilson,  and  Bieging  1983  (resolution  ~  T).  These  Te  estimates 
are  thought  to  be  less  sensitive  to  non-LTE  effects.  In  any  case, 
they  do  not  differ  from  those  determined  from  the  HI  10a  line 
observations  by  more  than  ~  30%  which  is  in  agreement  with 
the  finding  that  LTE  electron  temperatures  are  frequently  close 
to  actual  Te  (see,  e.g.,  Wilson  1980). 

Sources  B-E  have  line  velocities  which  are  all  within  6.5  km 
s-1  of  each  other.  Based  on  complimentary  H2CO  absorption 
measurements  (and  earlier  H  I  observations),  Downes  et  al. 


TABLE  1 

Recombination  Line  Measurements  near  S53 
A.  Results  from  HI  10a  Observations  (Downes  et  al.  1980) 


Continuum  Line  Emission  Electron 


Temperature*  Temperature*  Measure  Temperatureb 

Tc  Tl  Velocity  EM  Te 

Source  R.A.  (1950)  Decl.  (1950)  (K)  (K)  (kms-1)  (cm-6  pc)  (K) 


A .  lg^"^*  — 13°12'00"  3.2  0.15  32.3  3.2  x10s  7300 

B .  18  22  41  -13  18  31  3.2  0.34  45.5  4.4  x  104  6200 

C .  18  22  12  -13  17  41  9.5  0.68  52.0  5.9  x  10s  5800 

D .  18  22  28  -13  11  50  4.0  0.20  47.5  4.4  x  104  5300 

E .  18  22  15  -13  12  02  3.2  0.17  47.0  2.9  x  104  7600 


B.  Results  from  H76a  Observations  (Wink  et  al.  1983) 


Source 

R.A. (1950) 

Decl.  (1950) 

Line-to-Continuum 

Ratio 

TJTC 

Velocity 
(km  s-1) 

Electron 

Temperature 

Te 

(K) 

A . 

18h22m535 

— 13°12'00" 

0.232  +  0.019 

33.7  +  1.3 

5100  +  600 

C . 

18  22  12 

-13  1741 

0.240  ±  0.022 

53.0  ±  1.3 

6300  +  700 

*  Brightness  temperature  =  antenna  temperature/0.65. 

b  The  errors  in  the  electron  temperatures  are  estimated  by  Downes  et  al.  1980  to  be  +  2500  K. 
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TABLE  2 


VLA  Results  for  Sources  Near  S53 


90  Centimeters  20  Centimeters 


Integrated  Integrated  Spectral 

Peak  Flux  Density  Peak  Flux  Density  lNDEXb 

Brightness  S90cm  Brightness  S20cm  a 

Source  R.A.  (1950)“  Decl.  (1950)*  (Jy  beam  *)  (Jy)  (Jy  beam  l)  (Jy)  (90-20  cm)  Size* 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9) 


A .  18h22"53!  — 13°12'00"  <0.05  ...  0.97  1.07  +  0.05  >+2.0  <0!7 

B .  18  22  37  -13  19  00  0.48  0.61  +0.20  0.50  3.20  +  0.40  +1.1  4.4  x  2.8 

C .  18  22  13  -13  1740  1.85  2.61  +0.20  2.22  6.30  +  0.80  +0.6  2.4  x  1.3 

D .  18  22  29  -13  11  30  1.40  1.84  +  0.06  0.53  4.07  +  0.90  +0.5  4.2  x  3.3 

E .  18  22  09  -13  1140  0.78  2.77  +  0.11  0.27  >1.80  >-0.3  7.5  x  3.3 

F .  18  21  52  -13  1049  2.00  5.69  +  0.27  0.44  1.24  +  0.28  -1.0  3.6  x  0.8 


*  Positions  are  obtained  from  the  20  cm  map. 
b  S  oc  v+“. 

c  Sizes  are  estimated  from  the  20  cm  map  (Fig.  1)  expect  for  source  E  which,  because  of  its  large  size  and  low  surface  brightness,  was  best  estimated 
from  the  90  cm  map. 


(1980)  conclude  that  sources  B-E  are  all  at  the  near  kinematic 
distance  of  ~4.5  kpc.  We  therefore  accept  these  sources  as 
being  physically  associated  with  the  S53  complex  since  the 
brightest  source  C  is  coincident  with  the  H  ii  region  Sharpless 
S53  (Sharpless  1959).  However,  the  analysis  which  follows  is 
independent  of  any  physical  association  between  the  sources. 
Source  A  also  appears  to  be  an  H  n  region  but,  because  of  its 
large  velocity  difference,  is  presumably  not  physically  associ¬ 
ated  with  the  S53  complex. 

Table  2  summarizes  the  results  from  our  VLA  maps.  Source 
sizes  were  obtained  from  the  higher  resolution  20  cm  map 
except  for  source  E  whose  large  size  was  better  estimated  from 
the  90  cm  map.  Comparison  of  Figure  1  with  Figure  2  shows 
that  source  A  is  undetected  at  90  cm,  presumably  because  of 
self-absorption. 


Source  F  is  stronger  at  90  cm  than  at  20  cm  and  has  a 
nonthermal  spectrum.  Odegard  (1986)  has  suggested,  on  the 
basis  of  57.5  MHz  observations  made  at  the  Clark  Lake  Radio 
Observatory  (Erickson,  Mahoney,  and  Erb  1982),  that  it  is  a 
previously  unidentified  Galactic  supernova  remnant 
(G18.1  —0.2).  It  could  thus  be  associated  with  S53  but  since  we 
are  only  concerned  here  with  determining  H  n  region  proper¬ 
ties,  we  do  not  discuss  it  further. 

Integrated  flux  densities  are  also  listed  in  Table  2.  These 
were  obtained  by  averaging  a  fitted  Gaussian  model  value  and 
a  direct  contour  map  integration  value  with  errors  conserva¬ 
tively  defined  to  be  the  difference  between  the  two  determi¬ 
nations.  In  all  cases  these  error  estimates  exceed  those  expected 
due  to  random  noise. 

Figure  3  shows  an  overlay  of  our  20  cm  VLA  map  (Fig.  1) 


I8h  23m00s  22m30s  22mOOs  2lm30s 

RIGHT  ASCENSION 


Fig.  3. — Overlay  of  the  6  cm  single-dish  continuum  map  of  AltenhofT  et  al.  (1978)  with  our  20  cm  VLA  map  shown  in  Fig.  1.  The  half-power  beam  widths  for  the 
VLA  and  the  Effelsberg  100  m  telescope  are  illustrated  in  the  lower  right-hand  corner. 
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TABLE  3 

Flux  Densities  for  H  ii  Regions  Near  S53  from  Centimeter  Wave¬ 
length  Continuum  Observations 


Source 

Wavelength 

(cm) 

Peak 

Brightness 
(Jy  beam1) 

Integrated 
Flux  Density 

S 

(Jy) 

References 

A . 

11 

0.93  +  0.56 

1 

6 

1.35 

1.50  +  0.18 

2 

3 

1.40  +  0.36* 

3 

B . 

11 

1.86  +  0.56 

1 

6 

1.36 

2.90  ±  0.30 

2 

C . 

11 

5.28  +  0.56 

1 

6 

4.70 

5.50  +  0.18 

2 

3 

4.22 

4.75  +  0.36 

3 

D . 

11 

3.35  +  0.56 

1 

6 

1.70 

5.10  +  0.34 

2 

E . 

6 

1.60 

8.80  +  0.60 

2 

*  Estimated  from  the  Handa  et  al.  1988  map  assuming  that  source  A  is 
unresolved  by  their  observing  beam. 

References. — (1)  Reich  et  al.  1984;  (2)  Altenhoff  et  al.  1978;  (3)  Handa  et  al. 
1988. 

with  the  6  cm,  single-dish  continuum  map  of  Altenhoff  et  al. 
(1978).  The  Altenhoff  et  al.  map  was  made  with  a  beamwidth  of 
~2(7.  As  expected,  all  of  the  flat  spectrum,  thermal  sources 
A-E  clearly  appear  on  the  6  cm  map  while  the  steeper  spec¬ 
trum,  nonthermal  source  F  becomes  less  prominent.  The  1 1  cm 
(resolution  ~43;  Reich  et  al.  1984)  and  3  cm  (resolution  —  310; 
Handa  et  al.  1988)  single-dish  Galactic  plane  surveys  show 
very  similar  structure  for  the  area.  Where  possible,  integrated 
flux  densities  for  each  of  the  five  thermal  sources  A-E  have 
been  collected  in  Table  3  from  these  single-dish  surveys.  These 
values  are  taken  from  the  published  source  lists  in  the  original 
references  so  that  not  all  sources  have  listed  flux  densities  at  all 
frequencies. 

Errors  for  the  integrated  flux  densities  listed  in  Table  3  are 
difficult  to  obtain  since  estimates  are  generally  not  given  in  the 
published  source  lists.  Therefore,  we  assume  that  the  principal 
source  of  error  arises  in  the  difficulty  of  estimating  the  local 
background  levels  on  the  survey  maps  and,  somewhat  arbi¬ 
trarily,  we  assign  an  error  of  one-half  of  the  assumed  local 
background  brightness.  This  is  very  conservative  since  back¬ 
ground  levels  have  probably  been  determined  more  accurately 
than  this  in  each  of  the  three  surveys. 

IV.  DISCUSSION 

The  radio  flux  density  of  a  spherical  H  n  region  in  the  LTE, 
homogeneous  nebula  approximation  (Hjellming,  Andrews,  and 
Sejnowski  1969)  is  given  by 

S  =  3.07  x  lO~2Tev2Q(l  —  e~r) ,  (1) 

where  S  is  the  integrated  flux  density  in  Jy,  Te  is  the  electron 
temperature  in  Kelvin,  v  is  the  observing  frequency  in  MHz,  Q 
is  the  solid  angle  subtended  by  the  source  in  steradians,  and  z  is 
the  optical  depth.  At  radio  wavelengths  z  is  given  (Mezger  and 
Henderson  1967)  approximately  by 

t  ~  1.643  x  105v~21  EM  Te~135  ,  (2) 

where  EM  is  the  emission  measure  in  cm~6  pc,  defined  as 


with  ne  the  thermal  electron  density  in  cm-3  and  l  the  path 
length  in  parsecs  along  a  line  of  sight  L  through  the  source. 

In  the  optically  thin  limit  of  very  low  opacity  (z  1 ; 
1  —  e~T  ~  t)  equations  (1)  and  (2)  reduce  to 

S  =  3.07  X  lO~2Tev2nz  =  5  X  103re-o-35v“o  ln  EM  (4) 

and  in  the  optically  thick  limit  of  very  high  opacity  (t  1 ; 
1  _  e~z  ~  1)  equations  (1)  and  (2)  give 

S  =  3.07  x  10“27;v2fi  .  (5) 

The  short-wavelength-determined  Te,  EM,  and  S  values  pre¬ 
sented  in  Tables  1  and  3  combined  with  equations  (1)  and  (2), 
in  principle,  entirely  determine  the  spectra  for  the  five  H  ii 
regions  A-E  if  the  sources  are  uniform  and  the  filling  factors 
are  known.  In  practice,  effective  source  sizes  (filling  factors)  are 
poorly  determined  and  the  EM  obtained  from  recombination 
line  measurements  are  only  lower  limits.  However,  since  our 
90  cm  measurements  (Table  2)  are  in  the  optically  thick 
portion  of  the  spectrum  of  the  H  n  regions,  we  can  now  con¬ 
strain  the  source  parameters  far  better  than  previously  pos¬ 
sible. 

First,  let  us  consider  the  previously  available,  short- 
wavelength-determined  information.  Assuming  that  all 
sources  are  optically  thin  at  11, 6,  and  3  cm  (an  approximation 
which  is  justified  by  estimates  of  t  from  eq.  [2]),  equation  (4) 
can  be  solved  for  the  effective  source  size  as 

D  =  5  x  10-%cmT/-35  EM"1  ,  (6) 

where  S3  cm  is  the  3  cm  flux  density  in  Jy.  Our  20  cm  VLA  flux 
density  values  are  not  used  since  that  wavelength  lies  in  an 
intermediate  range  where  the  z  1  approximation  may  not  be 
valid.  To  obtain  an  accurate  S3  cm  and  minimize  the  problems 
in  determining  baseline  levels  inherent  in  single-dish  Galactic 
surveys,  we  fit  a  S  oc  v~01  spectrum  to  all  11,  6,  and  3  cm 
measurements  available  (Table  3)  for  each  H  ii  region,  weigh¬ 
ting  all  wavelengths  equally.  These  “best”  values  for  S3  cm 
obtained  from  the  fit  (Table  4,  col.  [2]),  together  with  Te  and 
EM  from  Table  1  (repeated  in  Table  4,  Cols.  [3]  and  [4]  for 
convenience),  were  then  used  in  equation  (6)  to  calculate  short 
wavelength  effective  source  sizes  QSA  (Table  4,  col.  [5]).  From 
these  short-wavelength  results  alone,  an  entire  spectrum  for 
each  H  n  region,  including  the  optically  thick  part,  is  now 
determined.  Substituting  the  values  from  Table  4  into  equa¬ 
tions  (1)  and  (2)  yields  the  dashed  curves  shown  in  Figures 

TABLE  4 

Parameters  Determined  from  Short-Wavelength 
(Optically  Thin)  Radio  Measurements 

Effective 

S3  om“  7'eb  EMb  Source  Size 


Source  (Jy)  (K)  (cm  6  pc)  Slsi  (sr) 
(1)  (2)  (3)  (4)  (5) 


A .  1.21  5100  3.2  x  105  3.7  x  10“8 

B .  2.17  6200  4.4  x  104  5.2  x  10“7 

C .  4.84  6300  5.9  x  105  8.7  x  10“ 8 

D .  3.85  5300  4.4  x  104  8.8  x  10“ 7 

E .  8.21  7600  2.9  x  104  3.2  x  10“6 


*  Estimated  by  fitting  an  a  =  —0.1  spectrum  to  the  11,  6, 
and  3  cm  values  (when  available)  from  short-wavelength 
surveys. 

b  Obtained  from  recombination  line  measurements 
(Table  1). 


©  American  Astronomical  Society  •  Provided  by  the  NASA  Astrophysics  Data  System 


I08  I09  io10 

FREQUENCY  1/  (Hz) 


Fig.  4 a 


FREQUENCY  v  (Hz) 
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Fig.  4. — Radio  spectra  for  the  five  H  H  regions  A-E  (see  Fig.  1)  near  the  S53  complex.  Data  points  are  taken  from  Tables  2  and  3  with  the  best  estimates  for  their 
appropriate  errors  shown.  The  dashed  curves  are  the  spectra  expected  if  the  source  parameters  are  entirely  determined  by  existing  short  wavelength  observations. 
The  solid  curves  are  those  obtained  if  the  parameters  of  the  FI  ii  regions  are  refined  through  use  of  our  90  cm  VLA  observations.  As  is  clear,  the  calculated  spectrum’s 
fit  to  the  data,  and  therefore  the  physical  description  of  the  regions,  is  greatly  improved  by  our  new  data. 
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4a-4e.  On  these  figures  we  have  also  plotted  all  flux  densities 
available  between  90  cm  and  3  cm  for  each  of  the  H  ii  regions 
(from  Tables  2  and  3)  including  estimates  (see  preceding 
section)  of  their  errors. 

It  is  clear  from  Figures  4u— 4e  that,  except  for  source  C, 
existing  (short-wavelength-determined)  parameters  predict 
continuum  flux  densities  at  90  cm  which  are  at  least  an  order 
of  magnitude  larger  than  those  actually  observed.  Such  dis¬ 
crepancies  can  be  explained  in  two  ways:  (1)  our  90  cm  flux 
densities  are  substantially  underestimated,  or  (2)  the  model 
parameters  determined  from  only  the  short-wavelength  mea¬ 
surements  are  significantly  in  error.  The  former  is  unlikely  (see 
§  II),  so  let  us  consider  the  latter. 

Any  error  in  the  models  must  lie  in  the  parameters  fi,  Te, 
and/or  EM  determined  from  the  short  wavelength  measure¬ 
ments  alone,  or  in  the  value  of  S3cm.  Also,  any  error  must  be 
very  large  since  the  predicted  (Figs.  4a-4e,  dashed  lines)  and 
measured  90  cm  flux  densities  generally  differ  by  more  than  an 
order  of  magnitude. 

The  value  of  S3  cm  is  unlikely  to  be  in  error  by  as  much  as 
50%  and  is  almost  certainly  not  the  origin  of  the  order-of- 
magnitude  discrepancies  we  find.  Errors  in  its  measurement 
come  mainly  from  uncertainties  in  determining  background 
levels  associated  with  the  extended  Galactic  nonthermal  emis¬ 
sion.  However,  our  “best”  estimates  for  S3cm,  obtained  from 
the  fitting  procedure  described  above,  are  quite  consistent  with 
slightly  longer  wavelength  values  and  with  direct  measure¬ 
ments  (when  available)  at  3  cm  (Handa  et  al.  1988). 

Values  for  Te  obtained  from  the  HI  10a  recombination  line 
observations  (Table  1)  assume  LTE  but  have  typical  quoted 
errors  of  only  ~  10%-40%.  Moreover,  results  of  the  H76a  line 
survey  (which  is  relatively  insensitive  to  non-LTE  effects) 
available  for  sources  A  and  C  (see  Table  1)  compared  with 
those  of  the  HI  10a  survey  show  that  departures  from  LTE 
electron  temperatures  are  not  more  than  ~30%.  That  LTE 


values  for  Te  are  approximately  correct  has  more  generally 
been  found  to  be  the  case  (see,  e.g.,  Wilson  1980)  and  is  to  be 
expected  since  the  non-LTE  effects  of  stimulated  emission  and 
of  pressure  broadening  tend  to  cancel  each  other  (Seaton  1980). 
Also,  non-LTE  effects  tend  to  underestimate  the  true  Te  which 
would  increase  the  discrepancy  between  the  measured  and  pre¬ 
dicted  values  of  Sgo  cm  rather  than  decrease  it. 

It  appears,  therefore,  that  the  large  differences  between 
model  predictions  and  measurements  must  arise  in  the  short- 
wavelength  estimates  of  EM  and  their  interplay  with  effective 
source  sizes  (eq.  [6]).  This  is  not  unexpected  since  estimates  of 
EM  from  recombination  line  observations  are  known  to  be 
lower  limits. 

Accepting  Te,  Sgo  cm,  and  S3  cm  as  relatively  well  determined 
and  sources  A-E  as  optically  thin  (t  <i  1)  at  3  cm  and  optically 
thick  (t  >  1)  at  90  cm,  by  assuming  that  the  effective  source 
size  has  no  frequency  dependence  (fi90  cm  =  fi3  cm)  (the  validity 
and  effects  of  this  assumption  are  discussed  below),  we  can 
solve  equations  (4)  and  (5)  to  give 

EM  =  l-65Tf,1'35S3cmS90cm_1  .  (7) 

These  new  estimates  for  the  EM  are  listed  in  Table  5,  column 

(2). 

Figures  4a-4e  show  the  spectra  calculated  from  models 
where  the  EM  values  from  Table  5  have  been  used  ( solid 
curves).  As  they  must,  these  curves  pass  through  the  “best” 
S3  cm  and  the  measured  S90  cm  values.  They  also  describe  well 
our  intermediate  frequency  S20cm  values  and  the  short- 
wavelength,  optically  thin  flux  density  points  from  the  liter¬ 
ature.  It  should  be  noted  that  these  new  EM  (Table  5,  col  [2]) 
are  significantly  different  from  the  previously  available  esti¬ 
mates  (Table  4,  col.  [4]). 

Physically,  determining  EM  based  on  Te,  Sgo  cin,  and  S3  cm 
(eq.  [7])  means  that  we  have  determined  an  effective  source  size 
fi  from  Sgo  cm  and  Te  rather  than  from  measurements  in  the 
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TABLE  5 

New  Parameter  Estimates  for  Sources  Near  S53 


EM*  £290crab  Qmap°  Filling  Factor 

Source  (cm  6  pc)  (sr)  (sr)  fi90 


A .  4.2  x  106  <3.0  xlO"9  <4.6  x  10  8  0.07 

B .  8.0  x10s  3.0  x  10~8  1.1  x  10~6  0.03 

C .  4.2  x  10s  1.3  x  10“7  2.6  x  10~7  0.50 

D .  3.8  x  105  1.1  x  KT7  1.2  x  10“6  0.09 

E .  8.7  x  10s  1.1  x  10-7  2.1  x  10-6  0.05 


*  Calculated  with  eq.  (7)  by  using  our  new  90  cm  integrated  flux  density 
measurements  (S90 cm ;  Table  2,  col.  [5]),  our  “best”  estimated  3  cm  flux 
density  measurements  (S3cm;  Table  4,  col.  [2]),  recombination  line  derived 
electron  temperatures  (Te;  Table  4,  col.  [3]),  and  by  assuming  that  the  sources 
are  optically  thin  (r  <  1)  at  3  cm,  optically  thick  (r  >  1)  at  90  cm,  and  have 
effective  source  sizes  which  are  frequency-independent  (fi3  cm  =  O90  cm)- 
b  Calculated  with  eq.  (5)  using  S90  cm  and  recombination-line-derived  Te. 
c  Calculated  from  measured  sizes  given  in  Table  2,  col.  (9);  see  note  c  to 
Table  2. 

optically  thin  portion  of  the  radio  spectrum  as  we  did  pre¬ 
viously  with  equation  (6).  It  is  useful,  therefore,  to  calculate 
these  new  values  of  fi90  cm  (Table  5,  col.  [3])  and  compare  them 
with  our  previous  short  wavelength  determined  estimates  (fiSA ; 
Table  4,  col.  [5])  and  with  direct  measurements  of  apparent 
source  sizes  on  our  20  cm  map,  Figure  1  (Hmap;  Table  5,  col. 
[4])- 

First,  it  is  apparent  that  the  measured  20  cm  sizes  (fimap)  are 
the  largest,  or  almost  the  largest,  in  all  cases.  This  is  to  be 
expected  since  the  fi  in  our  equations  is  not  strictly  the  size 
subtended  by  the  source  but  is  an  effective  size  which  includes  a 
filling  factor.  Second,  the  source  sizes  calculated  from  short- 
wavelength  measurements  are  inversely  proportional  to 
the  estimated  short  wavelength  EM  which  we  have  shown  to 
be  grossly  underestimated.  This  leads  to  size  estimates  which 
are  consistently  too  large  except  for  source  C  where  the  short- 
wavelength  estimate  appears  to  already  give  a  reasonable 
value. 

The  best  estimates  for  effective  source  sizes  are  fi90  cm  from 
our  observations.  These  can  be  used  to  obtain  source  filling 
factors  by  taking  their  ratio  to  the  apparent  physical  source 
sizes  Qmap.  This  ratio  is  given  in  Table  5,  column  (5).  (Note  that 
these  filling  factors  are  only  limits  in  cases  where  the  source  is 
unresolved  at  20  cm;  i.e.,  for  source  A.) 

The  filling  factors  in  Table  5  range  from  0.03  to  0.50  and 
indicate  that  significant  clumping  is  present  in  all  sources.  How 
does  this  affect  our  basic  assumption  that  the  gas  which  we  see 
in  the  optically  thin  regime  is  the  same  as  that  which  we  see  in 
the  optically  thick  regime — i.e.,  that  fi90cm  =  fi3cm?  Let  us 
consider  the  worst  case  of  an  H  n  region  consisting  of  high- 
density  clumps  embedded  in  a  lower  density  extended  medium 
such  as  a  raisin  pudding  model. 

At  first  examination  such  a  scenario  would  appear  to  cause 
difficulties.  At  shorter  wavelengths  where  both  components  are 
optically  thin  most  of  the  emission  arises  from  the  gas  with  the 
highest  EM,  i.e.,  from  the  clumps,  while  at  the  longer  wave¬ 


lengths  the  clumps  become  optically  thick  first  and  emission 
from  the  diffuse  component  dominates.  Such  considerations 
imply  that  source  sizes  determined  from  long  wavelength  mea¬ 
surements  should  tend  to  overestimate  Si  for  the  clumps. 
However,  comparison  of  Slgo  cm  (Table  5,  col.  [3])  with 
(Table  4,  col.  [5])  shows  that  to  clearly  not  be  the  case.  This 
also  implies  that  the  diffuse  gas  (the  pudding),  if  present,  is  still 
optically  thin  at  90  cm. 

Although  the  possibility  exists  that  our  new  estimates  are 
only  lower  limits  to  the  EM  and  upper  limits  to  the  effective 
source  sizes,  they  appear  to  be  the  best  values  currently  avail¬ 
able  for  H  II  regions  near  the  S53  complex. 

V.  CONCLUSIONS 

We  have  used  20  and  90  cm  VLA  observations  of  five  H  ii 
regions,  four  of  which  may  be  physically  associated  with  the 
S53  complex,  to  measure  their  continuum  flux  densities  and 
apparent  angular  sizes.  We  have  then  combined  these  data 
with  electron  temperatures  and  flux  densities  determined  from 
centimeter  wavelength  recombination  line  and  continuum 
observations  to  obtain  improved  estimates  of  the  emission 
measures  and  filling  factors  for  each  of  the  H  ii  regions. 

Previously  available  emission  measure  estimates  have  been 
shown  to  often  be  gross  underestimates  by  an  order  of  magni¬ 
tude  or  more.  Effective  source  sizes  determined  by  our  90  cm 
observations  compared  with  apparent  source  sizes  on  our 
20  cm  map  show  that  the  filling  factors  for  these  H  ii  regions 
range  from  -  0.03  to  0.5  and  indicate  that  the  gas  is  signifi¬ 
cantly  clumped. 

Our  new  results,  combined  with  source  modeling,  provide 
greatly  improved  estimates  of  H  n  region  physical  parameters 
and  allow  accurate  prediction  of  their  radio  spectra  over  more 
than  15  octaves. 

VI.  SUGGESTIONS  FOR  FUTURE  WORK 

The  question  of  clumping  in  H  ii  regions  is  important  since 
it  affects  all  parameter  estimates.  To  investigate  this  further,  we 
recommend  that  a  search  be  made  for  low-frequency  recombi¬ 
nation  lines  from  the  S53  complex.  These  lines  are  produced  by 
stimulated  emission  and  cannot  be  emitted  from  high-density 
gas  because  of  opacity  and  pressure  broadening.  Therefore, 
they  become  an  ideal  probe  for  the  presence  of  any  low- 
density,  distributed  gas  (the  pudding  around  the  raisins).  Also, 
ideal  conditions  are  present  for  production  of  these  lines  in  the 
direction  of  S53  because  the  presence  of  strong,  nonthermal 
background  emission  from  the  inner  Galactic  plane  can 
enhance  stimulated  emission. 

That  low-frequency  recombination  line  observations  should 
enable  the  construction  of  more  detailed  models  of  clumpy  H  ii 
complexes  is  not  a  new  idea  (see  Anantharamaiah  1986  and 
references  therein),  but  little  work  has  so  far  been  done.  Partly 
this  is  due  to  the  fact  that  instruments  with  sufficient  resolution 
for  low-frequency  spectral  line  studies  of  individual  sources 
have  only  recently  become  available. 
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